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Abstract 

Highly dendritic graphene crystals up to 0.25 mm in diameter are synthesized by low pres- 
sure chemical vapor deposition inside a copper enclosure. With their sixfold symmetry and 
fractal-like shape, the crystals resemble snowflakes. The evolution of the dendritic growth 
features is investigated for different growth conditions and surface diffusion is found to be the 
growth-limiting step responsible for the formation of dendrites. The electronic properties of 
the dendritic crystals are examined down to sub-Kelvin temperatures, showing a mobility of 
up to 6300 cm 2 V _1 s _1 and quantum Hall oscillations are observed above 4T. These results 
demonstrate the high quality of the transport properties despite their rough dendritic edges. 

The advent of large-scale graphene grown by chemical vapor deposition (CVD) on transi- 
tion metals opens a viable and promising route towards the commercialization of graphene-based 
electronics. The growth of graphene on copper has attracted considerable interest due to the 
simplicity, scalability, affordability, and homogeneity of the synthesized film. While this method 
solves the obvious problem of small-scale production associated with exfoliated graphene, it of- 
ten results in a graphene film with lower electronic performance.^ Significant efforts have been 
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made to reduce the extrinsic performance-limiting factors such as chemical impurities and struc- 
tural damages, recently, several studies have focused on improving the intrinsic electrical 
properties of CVD graphene. In particular, theoretical and experimental works have identified 
grain boundaries as one of the main sources of disorder in CVD graphene films. ^HEl 

Two approaches have been considered to overcome this problem. One consists in improving 
the electronic transport through the grain boundaries by engineering the growth conditions.^ The 
second strategy aims at decreasing the number of nucleation sites and increasing domain size in 
order to reduce the impact of grain boundaries on the electrical properties of the film. Follow- 
ing the pioneering work by Li et al.J^- 4 several CVD processes have been proposed^^HIIl i grow 
large crystals with lateral lengths up to 2.3 nimP These crystals display various morphologies: 
hexagons, flowers, squares and dendritic hexagons. Here, we report the growth of large, highly 
dendritic graphene crystals which we dubbed graphlocons due to their resemblance to snowflakes. 
Monolayer graphlocons, up to 250 fim in lateral size with very few defects were grown. We 
compared their growth shape evolution with other large island growths in order to confirm their 
unique morphology and propose a mechanism for the formation of dendrites. Field-effect tran- 
sistors (FETs) were fabricated on Si02/Si based on graphlocons and field-effect mobilities up to 
6300 cm 2 /V _1 s _1 were measured at 4 K. These devices also displayed well-developed quantum 
Hall effect (QHE) features despite their dendritic edges. 

Results and discussion 

Graphlocons were synthesized inside a copper-foil enclosure by employing a technique similar to 
the one reported by Li et al.J^but using a vertical quartz tube and higher gas pressure. The Cu 
enclosure was first annealed at 1025°C for 30 min in 150 mTorr of H2 flowing at 3 seem. The 
growth was performed for 30 mins at 1025°C at a pressure of 1500 mTorr, using a 0.5 seem CH4 
flow and a 3 seem H2 flow. Once the growth is completed, graphlocons are visible optically by 
heating up the copper foil in air on a hot plate for about 2 min at 200°C. This simple procedure 
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results in the oxidation of the copper regions which are not covered with graphene, creating a 
high contrast with those protected by graphene.^ Using an optical microscope, most graphlocons 
appear as bright six-fold snowflakes over the colored poly crystalline copper substrate ([T^). The 
diameter of these domains varies between a few microns and 250 jim. Four-fold islands were also 
found in specific regions, indicating that the growth might be affected by the crystal orientation 
or morphology of the underlying copper substrate.^^ Moreover, we observed a much higher 
domain density in regions of rough copper surface. []J>, which shows the growth of graphlocons 
along a pre-existing scratch, demonstrates clearly that the copper morphology has an effect on the 
nucleation behavior.^ 




Figure 1: (a,b) Optical micrograph of as-grown graphlocons inside a Cu enclosure. The dark 
colored region corresponds to oxidized Cu.(c) Optical micrograph of a graphlocon transferred to a 
Si02/Si substrate and (d) Raman spectra taken on a branch (monolayer) and center (bilayer) of a 
graphlocon. The scale bars on (a,b,c) are 50 jim. 

The seeding role of copper impurities was also indicated by the presence of bilayer (and few- 



3 



layer) graphene at the center of some domains. These bilayer structures are easily observable once 
the domains are transferred onto a SiCVSi substrate. As [T]: indicates, more than one of those 
terraced structures could sometimes be seen in the central region, suggesting that graphlocons can 
be composed of more than one crystal. To confirm the presence of multilayer graphene, Raman 
spectra were taken on monolayer and multilayer regions of a graphlocon. The bottom spectrum 
of [Tjd was measured with the laser aiming in the middle of one of the branches. It corresponds to 
a graphene monolayer, with a 2D to G peak intensity ratio of hD^G ~ 2 and a 2D-peak FWHM 
of 35 cm -1 . The top spectrum was obtained by directing the laser on a darker central region. It 
yields a much smaller I2£>/Ig ratio (~0.7) and a 2D-peak twice as broad (FWHM = 63 cm -1 ), 
which indicates the presence of a bilayer/multilayer in the center of the graphlocon. 23 We also 
notice that in both spectra the defect-induced D-peak is very weak, indicating the high quality of 
the graphlocon. 

To highlight and quantify the distinct morphology of the graphlocons, we compared them to 
graphene islands resulting from other CVD methods employed for growing large crystals. Flower- 
shaped crystals were obtained using the vapor trapping method described by Zhang et al.^ and 
square-shaped islands were grown using conditions similar to those reported by Wang et al.^(see 
Supporting Information for the details of the growth). [2^ shows a log-log plot of the area (A) of 
individual islands versus their perimeter (P), as measured from SEM images. Data from all growth 
methods are included and the dashed lines associated with each type of growth are linear fits. The 
solid line corresponds to the relationship between perimeter and area for perfect hexagons. For all 
growth methods, islands grow with a scaling exponent a < 2 (A ~ P a ), as expected for branched 
or fractal growth like diffusion-limited aggregates. * 24 * 25 * The fact that a is lower for graphlocons 
(a = 1.43) than for other types of growths (a = 1.66) clearly demonstrate their higher degree of 
ramification. 

All growths presented in [2^ show a similar island shape evolution, which suggests that a single 
growth mechanism could be at work. A possible explanation was first proposed by Nie et al.^ 
who argued that graphene on Cu(l 1 1) is surface diffusion limited. In this growth regime, the shape 
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Figure 2: Shape evolution of graphlocons. (a,b) Plots of log (island area) as a function log (island 
perimeter) for various growth conditions and techniques. Each data point represents an individual 
island and SEM images are shown in (a) for some of them. Islands represented by red dots and 
black squares were grown using conditions similar to Ref. [16] and Ref. \TT\ respectively. The 
dotted lines are linear fits of data points with perimeter > 10 jim and the solid line corresponds to 
the behavior expected for perfect hexagonal (i.e. non-dendritic) islands, c) SEM image showing 
the effect of competitive capture on the graphlocon shape. The dotted line illustrate the capture 
zones and the arrows highlight the anisotropy of the growth. 
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of the island stems from two competitive processes: (1) carbon atoms or aggregates attach to the 
island boundary at a rate k — 1 /At, where At is the time between two atom impingements. (2) 
Atoms diffuse or detach/reattach along the island edge in order to preserve the thermodynamic 
shape of the graphene island. Assuming a random diffusion process, the time td needed for an 
atom to diffuse along a boundary of size L is td oc L 2 /£>, where D is a diffusion coefficient. When 
At > td, atoms have enough time to diffuse before a new impingement occurs and the equilibrium 
shape dominates. This shape, which minimizes the edge free energy, can be found by the Wulff 
constructionists! and corresponds to a compact hexagon with zigzag edges, as reported by several 
experiments.* 11 * 29 * 3 ^ As the island grows, the diffusion time increases and faster growing orienta- 
tions start "growing out" when td > At. For graphene on Cu(lll), the growth rate has a six-fold, 
"flower-like" symmetry, with slow and fast growing orientation corresponding to zigzag and arm- 
chair edges, respectively. * 26 * 27 * 3 ^ In this growth regime, dendrites arise from Mullins- Sekerka^ 
type shape instabilities and grow along faster growth orientations. 

A transition from compact to ramified morphologies can be seen in [2^ for graphene grown with 
the vapor trapping method. It corresponds to the point where the curve deviates from the hexago- 
nal geometry (a = 2) and starts following a dendritic growth (a < 2). This transition occurs when 
td = At, which defines the correlation length L c of the dendrites. A transition was not observed 
for the other types of growth due to the limited range of island sizes we synthesized. Interestingly, 
the extrapolated transition point of the graphlocon growth curve coincides with the vapor trapping 
one. Past this transition point, six-fold branches form and dendritic arms progressively grow on 
their sides. For large graphlocons, the development of secondary dendrites can be observed on 
primary dendrites, thus illustrating the self-similar nature of this growth. We note that dendrites 
grow preferentially with a 60° angle with respect to their parent dendrite (or branch), consistently 
with a six-fold growth symmetry. According to our interpretation of the growth, we should also 
expect a change in the shape evolution with the growth temperature since it affects the surface 
diffusion of carbon species (D can be described by an Arrhenius equation^). As indicates, 
growing graphlocons in the same conditions but at higher temperature changes the value of a from 
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1.43 (T = 1025°C) to 1.61 (T = 1050°C), which is consistent with an increase in D. Additional 
growth experiments show that dendrites can be suppressed at higher growth temperature (see Sup- 
porting Information). Furthermore, we observe that the island shape is affected by the proximity 
of neighboring islands such that branches tend to grow longer toward regions of low island density 
(|2J:). This is a clear hallmark of competitive capture between islands sharing the same diffusion 
field. ^EI] This competition for the capture of the same carbon species alters the capture zone of 
each island and results in an asymmetric growth rate. All aforementioned observations provide 
evidence that the growths investigated are surface diffusion limited. 

To assess the effect of dendrites on the electronic transport properties of graphlocons, we trans- 
ferred them onto a SiC>2 /Si substrate and electrically contacted their lobes to make a back-gated 
graphene FET. Two such devices were cooled down to 300 mK in a pumped 3 He refrigerator 
and magnetotransport measurements yielded similar results. Their morphological features are in- 
dicated in [2^. In what follows we only present data for one of them, shown in [3^. The sheet 
resistance was obtained by the Van der Pauw (VdP) methocP using the leads 1, 2, 3 and 5. 

[3J3 shows the change in sheet resistance (p) as a function of gate voltage (Vg) at 300 and 4 
K without magnetic field. During cool down, the position of the charge neutrality point (Voirac) 
shifted slightly from 7 V to 5.2 V and the overall resistivity decreases by 10%, indicating a metallic 
behavior. Both curves display an on/off ratio of ^9 within the gate voltage range displayed in [3^. 
To fit these curves we used a diffusive transport model similar to the one proposed by Morozov 
et alJ25l They showed that the inverse of the resistivity, after a contribution due to the short-range 
scattering ps is subtracted, depends linearly with gate voltage (p — Ps)~ l — I^C ox (Vg — Vb) + o reS9 
where ji is the field effect mobility, C ox the gate capacitance and o res the extrapolated residual 
conductivity at the charge neutrality point. However, in order to account for the difference in 
mobilities for holes and electrons as well as the existence of a residual density no at the charge 
neutrality point due to large scale inhomogeneities, we can write the total carrier density as n + p = 
\J (n — p) 2 + Hq and n — p — C ox (Vg — VDirac)/e, where n and p are the densities of electrons and 
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Figure 3: Electronic transport of a graphlocon FET on SiC^/Si. (a) Optical micrograph of the 
device. The scale bar is 30 jim. (b) Sheet resistance as a function of gate voltage at 300 K (red) 
and 4 K (blue), (c) Hall resistance and (d) Hall mobility as a function gate voltage at 0.3K. The 
dotted lines are fits to the data using a diffusive transport model. 
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holes, respectively.^ This leads to 

P = ^ \ \+Ps, (1) 

with C ox = 1 1.5 nF/cm. Using this equation, we extracted an electron mobility ji n ^4300 cm 2 V _1 s _1 
and ji p ^6300 cm 2 V _1 s _1 for holes at 4K. The residual density was found to be no ~ 3.9 x 10 11 
cm -2 and the short-range scattering resistivity ps ~105 £1. These values compare well with those 
commonly measured in exfoliated graphene. 

E6] For 

a more thorough investigation of the carrier mo- 
bility, we also measured the Hall resistance R# at low B-field tfSjp) and extracted the Hall mobility 
/jL H = Rn/p (|3]l). The gate voltage dependence of Rh for two carriers is given by 

Rr = (^-^) (2) 
(n-p) 2 + nl 

which agrees well with our measurements as shown in [3]:. The Hall mobility was found to vary 
significantly as a function of Vg, especially in the hole doped regime. The highest values of jin for 
holes and electrons match those of the field effect mobility /I defined above. 

The homogeneity and quality of the graphlocon sample is also reflected by the magnetotranport 
measurements which display clear quantum Hall physics. [4^ shows the longitudinal (R xx ) and Hall 
resistivity (R xy ) as a function of Vg measured in a perpendicular magnetic field B = 9 T and T = 
0.3 K in the sample shown in [3^. R xx and R xy were obtained by passing a small, low frequency 
current through contacts 3-5, and measuring the voltage between contact 1-2 and 2-4, respectively. 
R xx was multiplied by a geometrical factor of 4.5 as derived from the VdP method. The data shows 
clear quantum Hall features, with well-resolved Hall plateaus and deep minima of R xx at filling 
factors v = ±2, ±6 and ±10 as expected for monolayer graphene. In the inset of [4^, we show 
R xx as a function of gate voltage and magnetic field. The resulting Laudau fan diagram shows the 
emergence of the quantum Hall states. Quantized Landau levels appear as maxima lines coming 
out of B = T and their linear dependence in B and Vg agree with the behavior for monolayer 
graphene.^ 
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Figure 4: Magnetotransport of a graphlocon FET on SiC^/Si. (a) Hall resistance R xy (black line) 
and longitudinal resistance R xx (red line) as a function of gate voltage VGate at B = 9 T and T = 0.3 
K. The quantum Hall plateaus at v = 2, 6 and 10 are indicated. The inset shows the Laudau fan 
diagram of R xx as a function of gate voltage and perpendicular magnetic field, (b) R xx as a function 
of B measured at Voate = 80 V. The dotted line is a fit for the SdHO. The inset shows the Fourier 
transform of R xx vs B~ l and the peak represents the charge carrier density up to a factor 4e/h. 
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Finally, we investigated the Shubnikov-de Hass oscillation (SdHO) displayed by R xx at high 
gate voltage (Vg = 80 V) for B > 3.4 T, as shown in [4}}. The figure's inset shows the Fourier 
transform of R xx as a function of B _1 which displays a prominent peak at 62±5 T _1 . This value 
corresponds to the carrier density up to a factor 4e/h, yielding n = (6.0±0.5) x 10 12 cm -2 which is 
consistent with VG-V/); rac =79±6 V. Using the expression detailed by Babinski et al.,^the SdHO 
allows us to extract a quantum mobility jiq of 1100 cm 2 V _1 s _1 with the fit shown in^. The 
extracted quantum mobility is about five times smaller than the Hall and field effect mobilities and 
characterizes the effective broadening of the Landau levels due to disorder. Our observed magneto- 
transport features are comparable to those observed for typical exfoliated graphene samples^ and 
high-quality CVD-graphene samples. E2J40I shows that the fractal nature of graphlocons does 
not significantly alter the electronic properties of graphene, despite the importance of the edge in 
quantum Hall physics. 

Conclusion 

In summary, we have synthesized large, highly dendritic graphene islands named graphlocons by 
CVD inside a copper enclosure. By comparing the island shape evolution of graphlocons to other 
types of large island growths, we showed and quantified the distinct morphology of graphlocons. 
Based on this analysis we explained the formation of dendrites in CVD-grown graphene as the 
result of the competition between carbon attachment and diffusion along the graphene island in a 
surface-diffusion limited growth regime. Graphlocons were transferred onto Si02/Si, electrically 
contacted and a hole mobility as high as 6300 cm 2 V _1 s _1 was extracted. Similar Hall mobility 
values were found and magnetotransport measurements displayed well-developed QHE as well 
as strong SdHO. These observations all demonstrate the high quality of graphlocons and their 
potential for graphene-based electronics. 
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Supporting Information Available 
Other Growth Methods 

In addition to the growth method described in the article, we used two other growth techniques that 
are known for growing large graphene crystals by CVD on copper. The results of these growths are 
shown in[2^. Following the method described by Ref . [161 the vapor trapping growth was performed 
by placing a piece of copper foil inside a small quartz tube (|5^ and b). The foil was first annealed 
at 1050°C for 30 min in 50 mTorr of H2 flowing at 7 seem. Graphene islands were synthesized at 
1050°C at a pressure of 200 mTorr, using a 1 seem CH4 flow and a 12.5 seem H2 flow. 

The square-island growth (|5J: and d) was achieved by employing large H2/CH4 ratio as reported 
by Ref. [T71 The copper foil was first annealed during 30 minutes in a flow of 10 seem H2 at 1050°C. 
The growth was then performed at 1050°C in a flow of 50 seem H2 and 1 seem CH4 at a pressure 
of 770 mTorr. After the growth the sample was cooled rapidly to room temperature under H2 flow. 

High Temperature Growth 

Growths were performed at a temperature of ^1080°C. Copper foils were first pre-treated using 
acetone, acetic acid and isopropanol. The copper foil was then annealed for 65 minutes in a flow 
of 75 seem H2 at a pressure of 650 Torr. Following the annealing process the pressure was lowered 
to ^130 Torr and a flow of 0.15 seem CH4 was introduced for 30 minutes, while the H2 flow was 
maintained at 75 seem. After the growth the sample was removed from the oven and allowed to 
cool rapidly to room temperature under H2 flow. 

This material is available free of charge via the Internet at jhttp : //pubs . acs . org/ , 
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Figure 5: Other growth methods employed to grow large graphene islands, (a) Schematic of the 
vapor trapping setup with a copper foil inserted in a small inverted quartz tube and (b) SEM picture 
of a typical graphene island synthesized with this method, (c) Schematic of the square-island 
growth and (d) SEM picture of the synthesized islands. 
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Figure 6: Optical micrograph of graphene crystal grown using a high pressure, high temperature 
CVD process. The dark colored region corresponds to oxidized Cu 
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